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We present experimental and theoretical results of a detailed study of laser-induced continuum 
structures (LICS) in the photoionization continuum of helium out of the metastable state 2s 1 So- 
The continuum dressing with a 1064 nm laser, couples the same region of the continuum to the 
4s 1 So state. The experimental data, presented for a range of intensities, show pronounced ion- 
ization suppression (by as much as 70% with respect to the far-from-resonance value) as well as 
enhancement, in a Beutler-Fano resonance profile. This ionization suppression is a clear indica- 
tion of population trapping mediated by coupling to a contiuum. We present experimental results 
demonstrating the effect of pulse delay upon the LICS, and for the behavior of LICS for both weak 
and strong probe pulses. Simulations based upon numerical solution of the Schrodinger equation 
model the experimental results. The atomic parameters (Rabi frequencies and Stark shifts) are 
calculated using a simple model-potential method for the computation of the needed wavefunctions. 
The simulations of the LICS profiles are in excellent agreement with experiment. We also present 
an analytic formulation of pulsed LICS. We show that in the case of a probe pulse shorter than the 
dressing one the LICS profile is the convolution of the power spectra of the probe pulse with the 
usual Fano profile of stationary LICS. We discuss some consequences of deviation from steady-state 
theory. 
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I. INTRODUCTION 



A. Continuum coherence 



The properties of the continuum of quantum states observed in the photoionization of an atom or the dissociation 
of a molecule have attracted interest since the first formulations of quantum theory. Until the widespread use of 
lasers, the continuum was regarded as a dissipative system, to be treated by means of rate equations (typically using 
the approximation of the Fermi Golden Rule jlj to calculate the transition rate into the continuum). Contemporary 
views no longer regard continuum states as an irreversible drain of probability. An electron acted upon by a laser 
field of sufficiently high frequency to produce ionization does not move irreversibly away from the binding region 
|^| |J. Laser interactions with a continuum have been shown to exhibit Rabi-type oscillations (if the laser field is 
sufficiently strong) |,@] and to enable nearly complete population transfer into a continuum (i.e. photoionization or 
photodissociation |7|) or from a continuum (i.e. photorecombination |8)). Both theoretical j?],[3]-[l2| and experimental 
P, fn3|]i!| ] work has shown that the wavefunction of an electron in a laser field may often maintain coherence during 
the course of photoionization, and that photodissociation can also maintain coherence. Consequences of continuum 
coherence are to be found in the well-studied phenomena of autoionization p^j-pTj , which is revealed as a resonance in 
plots of photoionization cross-sections versus frequency. Such resonances display the effect of interference between two 
ionization channels whose destructive interference can lead to complete suppression of photoionization at a specific 
wavelength. 

One of the more interesting coherent phenomena involving the continuum is the use of a laser field to embed an 
autoionization resonance into an otherwise featureless photoionization continuum. This dressingof a continuum by a 
bound state by means of a laser field is termed a "laser induced continuum structure" (LICS) fll£fr-p5||. 

As with other types of resonances described by scattering theory, LICS can be observed in various reaction channels. 
In the LICS observed in photoionization a strong laser field, connecting an unpopulated bound state to the continuum, 
mixes some bound-state character into the continuum; it places a resonance into the continuum. This optically induced 
resonance can be detected by a second (possibly weaker) laser field which serves as a probe and induces a transition 
from a populated bound state into the newly structured continuum. 

We remark that when the continuum is composed of several channels (derived e.g. from different asymptotic internal 
states of dissociation products) the dissociation has been shown to give rise to a quantum control scheme (2(| in which 
the ratio between the asymptotic channels can be controlled by tuning the dissociating (or dressing) laser across a 
LICS resonance. 

Theoretical descriptions of continuum structure near autoionizing resonances have been presented p7|-|3l|] as have 
LICS in the absence of such resonances ]32| , [33| . Observations of LICS in photoionization near autoionization resonances 
p3| , pT|j34H3^ | as well as observations of LICS in photoionization to an otherwise unstructured continuum 1 18-25) have 
been reported recently. Under appropriate conditions this LICS can lead to the suppression of photoionization, by 
creating a coherent-superposition of bound states immune to photoionization (a trapped state, or dark state p6|). 
The occurrence of significant photoionization suppression signals the creation of a dark state; only with coherent 
excitation can this aspect of continuum coherence be revealed. 



B. The resonance profile 



When a coherent radiation field resonantly excites a two-state system, then a probing of either energy level by a 
weak probe field reveals a splitting (the Autler-Townes splitting) attributable to the dressing of bare atomic states 
by a strong field |57j . When the strong transition connects a bound (but unpopulated) state to a continuum, a weak 
probe will again reveal structure, the LICS. 

The theoretical basis for describing LICS has been presented in numerous articles (see for example | f38|]39| ), for 
cases in which the fields are constant-amplitude continuous-wave monochromatic, and for pulsed radiation. As such 
work has shown, the LICS of photoionization typically takes the form of an asymmetric window resonance, for which 
the frequency dependence of the photoionization cross-section a{u>) may be parameterized by the form of a typical 
scattering resonance in a background Jl7| , ^0t 

ctM = c{u) + (1) 

where x = {lo — oj#)/(F/2) is a dimensionless detuning of the light frequency lu away from the resonance value lur, 
measured in units of the resonance- width parameter (or loss rate) T. Equivalently, one may use a form akin to that 
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developed by Fano jlfj and Fano and Cooper jl6| to describe the Beutler-Fano profiles observed in autoionization or 
electron-scattering cross-sections, 



a(w) = CX6 M + & a 2 



1 



(2) 



where q is the dimensionless Fano q parameter. The terms C(u>) or <7b(u;) are slowly varying (or constant) contributions 
to the cross section. Appendix A describes simple estimators of the parameters when fitting experimental data to 
this functional form. 

We found it useful to use a variant of these parameterizations, 



cr{D) — er a Re 
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where D is the detuning of the probe laser, and and Tq are parameters discussed in section II C 



Although we have presented profile formulas appropriate to photoionization cross-sections, what is important for the 
present paper is the functional dependence of an experimentally observed quantity, here cr(w), upon a dimensionless 
detuning variable x. It often proves convenient to fit dimensionless observables with these parameterizations, in which 
case the parameters A, B and C, or a a and <Jb, become dimensionless. 



C. Present work 

In an earlier paper [5sl we presented experimental evidence for substantial suppression of photoionization by means 
of LICS. As we stressed, only with the use of very coherent laser pulses, together with an appropriate choice of 
atomic states, is it possible to observe pronounced suppression. Here we present further experimental results showing 
dependence of the LICS profile on various pulse properties. 

Our experimental results are for LICS induced in metastablc helium, by means of a dressing laser (1064 nm) that 
couples the 4s 1 So state into a continuum. The structure is observed by scanning the frequency of a probe laser which, 
at a wavelength of 294 nm, couples the 2s x So state into this same continuum. Figure |l| shows the relevant energy 
levels and transitions. 

We discuss numerical simulation of the experimental results, and show that with calculations of atomic parameters 
the observations are in excellent agreement with computations. We present an analytic description of pulsed LICS, 
emphasizing the difference between pulsed and steady excitation. These results reduce to expressions ([!]) and (^) 
under appropriate conditions. Our calculations of needed atomic parameters (polarizabilities, Stark shifts, etc.) are 
done using a simple effective-potential model of the active electron (see Appendix [b]). We describe this method, and 
we show that this model provides quantitative agreement with experiments. 
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II. EXPERIMENTAL 



Our observations of LICS in the photoionization of helium are carried out using an atomic beam of metastable 
helium atoms, which cross spatially overlapping beams of pulsed laser radiation. We control the duration of the 
pulses and the overlap in space and time. Figure |^ shows a schematic diagram of the experimental apparatus. 

A. The apparatus 

a. The atomic beam. A pulsed beam of metastable helium atoms is formed as helium, at a stagnation pressure 
of 1200 mbar, expands through a nozzle (General Valve, opening diameter 0.8 mm) and enters a region of pulsed 
discharge located 4 mm from the nozzle. A skimmer (diameter 0.8 mm) placed 40 mm from the nozzle collimates the 
atomic beam and separates the source chamber from the region of interaction and detection. 

b. The metastable atoms. The metastable states of helium used for our photoionization experiments are populated 
in a gas discharge, where the atoms are excited from the ground state Is by pulsed electron impact. To enhance 
efficiency and stability of the discharge, a filament, placed 10 mm downstream, provides the discharge with a continuous 
electron current of several mA. The peak current, produced by ions and electrons in the pulsed discharge, reaches 20 
mA. The electron source enhances the efficiency of the metastable production by at least two orders of magnitude 
and the breakdown voltage is reduced to 1/50 of the value needed for an unseeded discharge. 

c. Diagnostics of metastable population We deduce the population of metastable atoms by using a two-photon 
version of REMPI (resonant-enhanced multiphoton ionization) to produce photoions, with pulsed radiation from a 
broadband dye laser (LPD 3000, Lambda Physik). The ions are observed as an electric current. Figure || shows the 
REMPI transitions used to detect the 2s 1 So and 2s 3 Si metastable atoms. 

Figure [| shows the ionization signal of the metastable singlet and triplet states as a function of probe-laser wave- 
length. The REMPI transitions are strongly broadened by saturation. We also confirmed that the ionization is 
saturated. Populations are inferred from the heights of the ionization signal peaks. Typically 90% of the metastable 
population is found in the triplet state 2s 3 Si , and 10% is found in the singlet state 2s 1 So . 

d. The lasers. The pulsed atomic beam crosses the coinciding axes of two pulsed laser beams: a UV probe pulse 
and an IR dressing pulse. The probe laser field, coupling the metastable state 2s 1 So to the ionization continuum 
(see Fig. [l]), is derived from a single mode cw dye laser system operating at 587 nm. This radiation is amplified in a 
pulsed dye amplifier, pumped by the second harmonic of an injection-seeded Nd:YAG laser. The wavelength of the 
cw radiation is measured to an accuracy of AA/A = 2 x 10~ 6 in a Michelson-type wavemeter, using a He-Ne laser, 
stabilized on an iodine line, as a reference. The pulsed radiation is frequency-doubled in a BBO crystal, thereby 
providing (probe) laser pulses with a pulse width At = 2.3 ns (half- width at 1/e of I(t)) at A = 294 nm. 

The pulse of the dressing laser couples the state 4s 1 So to the ionization continuum and induces the continuum 
structure. This pulse, obtained from the fundamental frequency of the Nd:YAG laser (A = 1064 nm), has a pulse 
width At = 5.1 ns. 

A folded optical delay line is used to adjust the time delay between the dressing and probe pulses. After passing 
through the optical setup, pulse energies of up to 2 mJ for the probe laser and up to 250 mJ for the dressing laser are 
delivered to the atomic beam. The laser diameters at the atomic beam position are 0.5 mm for the probe and 3.5 mm 
for the dressing laser, yielding peak intensities of up to Ip ^ — 100 MW/cm 2 and if?' = 300 MW/cm 2 , respectively, 
(we use a superscript 1^ to denote the peak value of a time- varying pulse I(t).) Because of the large diameter of the 
dressing-laser beam, the variation of its intensity across the probe laser profile is small. 

e. Ionization detection. Ionized helium atoms are mass-selectively detected with a double-thickness micro-sphere 
plate (El Mul Technologies) after passing through a 30 cm time-of-fiight segment. The output current of the micro- 
sphere plate is amplified with fast broadband amplifiers and integrated in a boxcar gated integrator (EG & G 4121 
B). 

B. Results 

The continuum structure induced by the dressing laser is revealed in the ionization of the metastable state 2s 1 So 
produced by the probe laser. When the probe laser is tuned across the two-photon resonance between the states 2s 
1 So and 4s 1 So, we observe a strong and spectrally narrow feature in the ionization cross-section (see the following 
subsections). The structure is due to the laser-induced mixing of the bound state 4s with the ionization continuum. 
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1. Background signal 



At high laser intensities, two-photon ionization of the metastable triplet state (see Fig. forms a background, 
which has to be subtracted from the total ion signal. The background ionization was determined by tuning the UV 
laser away from the region where LICS is observed and measuring the ionization probability with and without the IR 
laser present. The difference of these two currents (the difference signal) is attributed to the two-photon (IR + UV) 
ionization of the triplet level, the population of which exceeds that of the singlet level by about an order of magnitude. 

Figure shows the measured difference signal, relative to the total ion signal, as a function of probe intensity. 
The probe laser frequency is here tuned also far away from the two photon resonance between states 2s and 
4s 1 So- The plotted signal is proportional to the ionization probability of the 2s 3 Si state. For comparison we show 
the calculated ionization probability for 2s 3 Si state caused by one photon at 294 nm and one photon at 1064 nm, 
using the theoretical ionization probability. Using wavefunctions described in Appendix [b] we obtained the result 
^triplet — 1-6 x 10~ 9 I p Id- Here and elsewhere laser intensities are expressed in W/cm 2 and the unit for the rate is 
s _1 . The calculated relative transition probability is in good agreement with the experimental data. 

From Fig. H we see that the background signal contributes only a few percent to the total ion sig nal for I p 0) < 10 
MW/cm 2 , while it adds more than 30% for the highest probe laser intensity used here. This background signal has 
been subtracted for all the data presented below. 



2. LICS for weak probe laser 

If the probe laser is weak, the shape of the continuum structure is dominated by the strong dressing laser, and the 
probe does not affect the internal structure of the atom. It is therefore possible to distinguish between the influence 
of the dressing and probing lasers. 

Figure ^| shows, for a weak probe laser, the observed structure in the ionization continuum of helium for various 
dressing-laser intensities (in each column) and for two pulse timings (in each row). The left-hand column of frames 
show the LICS produced by coincident laser pulses; the right-hand frames show LICS produced when the dressing 
pulse precedes the probe pulse (negative delay) by 5 ns. The dotted lines are results from numerical simulations, 



discussed in section III 



The continuum resonances exhibit typical asymmetric window-resonance structure as a function of frequency, such 
as are observed in autoionization. The LICS exhibit an enhancement of photoionization, taking largest value P max 
at frequency v maxi and a suppression of photoionization, a dip to P m in centered at frequency v m in- 

As can be seen, it is possible to control the photoionization probability by appropriately tuning the wavelength of 
the probe laser, thereby selecting the region of enhancement or suppression of ionization. Photoionization suppression 
indicates population trapping. We observe photoionization suppression as strong as a drop to 30% of the value far 
from resonance [p5| . 

The primary effects of an increase in dressing-laser power are a broadening and a shift of the LICS profile, both for 
coincident and delayed laser pulses. This is expected: the dressing laser intensity is responsible for the width-parameter 
r and the Stark shifts. 

For delayed laser pulses the LICS is not as strong as, and is spectrally narrower than, for coincident pulses. This 
is because when pulses are not coincident, the dressing laser is intensity is lower when the probe pulse samples the 
continuum structure: the LICS width is smaller and less pronounced when the dressing is weaker. 



3. LICS for strong probe laser 

For a strong probe laser, it is no longer possible to distinguish between the effect of the dressing and probing laser 
pulses. Both pulses contribute to the continuum structure. Whether the influence of the probe pulse is negligible or not 
can be estimated from the ionization probability of state 2s 1 So- Figure [?] shows the measured ionization probability of 
state 2s 1 So as a function of the probe laser intensity (in the absence of the dressing laser). The observed dependance 
is in good agreement with the theoretical prediction (Fig. 0, dashed curve), based on the calculated ionization rate 
T 2s , p (see Table 0). 

For probe intensities of several tens of MW/cm 2 the ionization is strongly saturated. In this regime the interaction 
induced by the probe laser cannot be considered as a small perturbation of the atomic structure. The observed profile 
will consequently deviate from the profiles of eqns. (|l|) and 
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Figure |8| compares the observed LICS for weak and strong probe pulses. The peak values of probe laser intensities 
are Ip — 0.5 MW/cm 2 (Fig. ||a) and Ip ^ = 40 MW/cm 2 (Fig. ||b). The peak intensity of the dressing laser was 
ll ' = 44 MW/cm 2 and 50 MW/cm 2 . An increase of the probe laser intensity reduces the peak value of the ionization, 
relative to the value far from resonance, from 140% to 110%. This reduction of the photoionization enhancement is 
the result of saturation: the atoms are rapidly photoionized for all frequencies except those where photoionization is 
suppressed. Because the strong probe pulse ionizes about 90% of the population of the metastable singlet state (see 
Fig. |7|), the maximum enhancement of the ionization probability can be no more than 10%. 



C. Fitting Experimental Data 



To parameterize our experimental results we first subtracted a linearly varying background fitted to the far-from- 
resonance values, and then performed a least-squares fit of the resulting data to a version of the Breit-Wigner or 
Fano-Cooper formula, Eqn. (|^). This parameterization introduces two contributions to the overall width parameter 
r = + Tq. The ionization rate of the dressing laser, r^, can, in principle, be computed using the traditional Fermi 
Golden Rule. The remaining term, Tq, represents an additional empirical rate; it expresses the effect of spontaneous 
emission, of fluctuations of intensity or phase, and other stochastic processes. Figure ^ shows an example of the 
experimental data (small squares) and the best analytic fit to this curve (dashed line). As can be seen, the data is 
very well described by the analytic expression. Table [j] presents values of the least squares parameters for undelayed 
pulses. 



TABLE I. Fits of the experiment, At = 
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Here Ip and l d are the peak intensities of the probe and dressing 
lasers in MW/cm 2 , q is the Fano q parameter, is the ioniza- 
tion rate from the dressing laser, T is the overall width parameter, 
r = T — r^, and Ttheory is the theoretical value for r^. All rates 
are in units of GHz. 

From several such fits we found the trend of the width parameters as a function of dressing intensity. Figure [l(] 
plots values of the parameter Tq and of the empirically determined r^. We show also the theoretical value r theory of 
the photoionization rate T^, calculated using the method described in section [b|. Theory and experiment are in very 
good agreement. 

Laser intensity fluctuations lead to fluctuations in the laser-induced Stark shifts. These fluctuations give an addi- 
tional contribution to the LICS width, 

j^Stark / qO q0 \ / a\ 

l - W4s,d _ ^s.d/ Flue (4) 



Here Sfj is the peak Stark shifts of state i induced by laser field j (see Section 

typical laser intensity fluctuations of A/l ^ = ± 10 % Id the additional broadening is TQ tark /2t: — 0.2 GHz. 



[Hi). With iy> = 75 MW/cm 2 and 
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III. NUMERICAL SIMULATION 



Usually laser induced continuum structures are deduced from the equations for the density matrix of an atom. This 
approach allows one to take into account coherence dephasing processes such as laser frequency fluctuations, collisions 
etc. But in the case of coherent laser pulses, with transform-limited bandwidths, one can significantly simplify the 
description of the process by using the time-dependent Schrodinger equation for the amplitudes of atomic states 
involved in the process. 



A. The main equations 



We consider an atom interacting with two laser fields: a probe laser with amplitude £ p (t), intensity I p {t) and carrier 
frequency lo Pi and a dressing laser with amplitude £d(t), intensity Id(t) and carrier frequency u)d- The total electric 
field is 



E(i) = Re Vj£j(t) CX P (—ib)jt) , 

j=P,d 



(5) 



where ej is a unit vector. The probe and dressing lasers are tuned close to the two-photon Raman resonance between 
states 1 and 2, and their carrier frequencies differ from exact resonance by the two-photon detuning D, given by 



fiD = fi{iu p -u Jd )-{E 2 -E 1 ). 



(6) 



We express the state vector |W (t)) of the atom as a combination of two particular discrete states 1 and 2, with 
energies E\ and E2, together with a sum over all other bound states and an integral over continuum states: 



|*(t))= £ C n {t)exp{-iE n t/h)\n) 

n=l,2 

+ A n (t)exp(-iE n t/K)\n) 

+ £ / dEA Ed {t)exp{-iEt/h)\E,j), 



(7) 
(8) 

(9) 



where \E,j) is the continuum state with energy E (quantum number j identifies different possible continua). 

As is customary P, p^|p7|j3^ , [4l| fhij we adiabatically eliminate all amplitudes except those of states 1 and 2. In so 
doing, we obtain effective interactions, and dynamic energy shifts, expressible in terms of a polarizability tensor p7[, 



QmnM = ^ 



(m|d|fc)(fc|d|n) 



Ek — E„ 



{m\d\E,j)(E,j\d\n) 



E — En 



huj 



(10) 



where d is the dipole-moment operator. The continuum integral includes the resonance energy E = E m + Twj, which 
is dealt with by breaking the integration into a nonsingular principal value part (denoted by V) and a resonant part 
(yielding the transition rate of a Fermi Golden Rule) , 



ReQ m „H =VJ2 



(m|d|fc)(fc|d|n) 

Ek — E n 



Too 



(m\d\E,j)(E,j\d\n) 



E - En 



flLO 



ImQ r . 



^n(m\d\E = E m + fkj,j){E = E m + hw,j\d\n) 



(11) 



The adiabatic elimination yields a two-state time-evolution equation, which may be written 



d 




1 


i— 

dt 


C 2 (t)_ 


~ 2 



2Si(t) - iTi(t) - iji -(i + q)fl(t) 

-(i + q)Q(t) 2S 2 (t) - iT 2 (t) - t l2 - 2D 



Ci(t) 
C 2 {t) 



(12) 



This equation, the basis for our simulation, requires dynamic Stark shifts Si{t), detuning D, photoionization rates 
Fj(t), spontaneous emission rates 74, an effective Rabi frequency O(t), defined below in Eqn. (|17|), and the Fano 
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parameter q. These various parameters can be obtained from atomic wavefunctions, as the following paragraphs 
show. 

We are interested in the total ionization probability after pulses cease, P. This can be determined if we know the 
populations of states 1 and 2 after the action of the laser pulses: 

P=l-|C 1 (oo)| 2 -|C 2 (oo)| 2 . (13) 

The LICS profile is taken as the variation of P with two-photon detuning D. 



B. Formulas for Atomic Parameters 



The theoretical ionization rates and Stark shifts are obtainable by summing partial rates: 

r,(t) = ]T rt 3 (i), s t (t) = s-W. ( 14 ) 

j.a ct 

where (i) is the ionization rate from state i to continuum j caused by probe (a = p) or dressing (a — d) lasers 

and hSi a (t) is the dynamic Stark shift of the energy of state i produced by laser a: 

Sfa(t) = K(t)\ 2 e„- [Q«(w a ) + Q«(-Wa)] <■ (16) 

The quantity 



n(*) = Ev r itW r w(*) ( 17 ) 

3 

is an effective Rabi frequency for the two-step transition from the state 1 to the state 2 with intermediate population 
of the continuum states j with energy E ~ Tilj p — E2 + E% ~ tkud — E\ + E-i. The Fano q parameter can be evaluated 
from the expression 
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qQ(t) = \£ P {t)B d {t)\ e p - [Q 12 {lj p ) + Q 12 {-cJ d )} <■ (18) 

From Eqns. (0) and (|J) it is apparent that q depends only on atomic properties. 

In principle, all of these quantities can be calculated if the atomic wavefunctions are known. 

We note that the Fano parameter q is basically the ratio of polarizability (such as occurs in a Raman transition) 
to the product of two dipole transition moments into the continuum. In systems with \q\ 3> 1 the Raman- type 
transitions dominate. The ionization probability is enhanced, when tuning the probe laser across the two-photon 
resonance, because Raman-type transitions open additional (multiphoton) ionization channels from the initial state 
to the ionization continuum. If q — no Raman-type transitions are present and no enhancement of ionization can 
be observed. 



C. Values for Atomic Parameters 



For the calculations we need the values of many transition matrix elements, both bound-bound and bound- 
continuum. Due to the relative simplicity of the two-electron helium atom there exist several methods, of varying 
complexity and accuracy, which can provide the needed wavefunctions |44|-f48||. However, these methods are usually 
optimized for the calculation of particular parameters (such as the ionization energy). 

Therefore, for our purposes it is desirable to use some simpler method. The simplest of these are based on the 
assumption that each electron moves in an effective central potential, and that therefore the atomic wavefunction is 
separable into a product (antisymmetrized) of individual electron orbitals, each of which has a spherical harmonic for 
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angular dependence and a radial dependence given by the solution to an ordinary differential equation (see Appendix 

We have found that a simple single-electron model potential satisfies requirements for speed and accuracy. Appendix 
[b| describes the mathematics of this method for computing matrix elements. As noted there, a single empirically 
determined parameter Xg fixes the potential, and hence all the wavefunctions, for each orbital angular momentum I. 

Table [0] presents the theoretical values we used in the simulation. 

TABLE II. Theoretical values for simulation 

q 0.73 

Tip 13.4 x I p Si p 

T 2p 1.9 x I p S 2p 

Tid Sid 

r 2d 73.7 x I d S 2d 



9.0 x Ip 
13.0 x I p 
70.0 x I d 
142.0 x I d 



Here I p and Id are the intensity of probe and dressing lasers expressed in W/cm 2 . The ionization rates and Stark 
shifts are expressed in s . These results are in good agreement with more elaborate computations carried out by 
Klystra, Paspalakis and Knight f49[| . The agreement with these more accurate results, and the successful fitting of 
our experimental data, support the validity of our atomic wavefunctions. 

The numerical simulation of the observed LICS-profilc includes experimentally determined intensity fluctuations 
and integration over the spatial probe laser intensity profile. These effects have to be taken into account, because 
position, bandwidth and shape of the LICS are strongly dependent on the laser power. 

To account for pulse-to-pulse variations of the intensity we average our computed ionization probabilities over a 
distribution of peak pulse intensities. We used a Monte-Carlo method to estimate the relevant integrals, assuming 
intensities evenly distributed between 0.9 I p and 1.1 I p . The results of the simulations are in excellent agreement with 
the experimental results, as are seen in Fig. o. 



D. Deviations from Steady State 



When the probe pulse is weak, and much shorter than the strong dressing pulse, the excitation can be regarded, 
to a first approximation, as an example of steady-state LICS (see sections IVD - IV F). One expects, under these 
circumstances, that an observed LICS profile will follow the form of Eqns. (|1[) or (|2|) 39j] or, equivalently (j^), and 
that the parameters qfu and r f it obtained from fitting such an expression will be equal to the values qtheo and r t ^ eo 
obtained by evaluating matrix elements. 

To show the effects of pulse shapes we display, in Fig. |ll| the profiles predicted in two approximations: 

• Steady-state profile, using the theoretical q value. 

• Simulation for a single pair of pulses (i.e. transient theory), again using the theoretical q value. 

The difference between these curves is small but it shows the importance of transient effects and the limitation 
of the steady-state approximation of sections IV D - IV F . Here "transient" means that we determine the ionization 
probability from the solution of the differential equations (|2|) thereby taking into account the pulsed variability of the 
probe interaction. The spectral width of this pulse, approximately l/r p , effectively contributes additional broadening, 
To, to the LICS profile. Only with the steady-state approximation (r p — > oo and Tq = 0) is there complete suppression 
of photoionization. 



E. Delayed Pulses 



In addition to the results just described, for which the probe pulse was timed to coincide with the peak value of the 
dressing pulse, we have also obtained results with probe pulses delayed with respect to the peak of the dressing pulse. 
As with the previously discussed results, it is possible to fit these LICS profiles with the function shown in Eqn. <^). 
Table III E gives the values of parameters from these fits. 
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TABLE III. Fits of experiment, At = -5 ns 
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98 



The width parameters obtained from fits for delayed pulses are significantly larger than what one would expect 
from the value of the dressing intensity at the peak of the probe pulse. To emphasize this point, Fig. [jj] shows the 
dressing and probe envelopes. The point marked B indicates the dressing intensity required to produce the observed 
width. The point marked A shows the (lower) intensity occurring at the peak of the probe pulse. As these points 
show, it is not possible to evaluate the LICS profile by considering only the moment of peak probe intensity. The 
observed profile averages in time the Stark shifts and the ionization rates. 

F. Summary 

In the simulations of the experimental data presented above, intensity fluctuations of both lasers as well as an 
integration over the spatial probe laser profile are taken into account to describe the observed LICS. In the present 
experiments the diameter of the dressing laser was much larger than that of the probe laser, and so the spatial 
variation of the dressing laser intensity was neglected in the simulation. 

a. Weak probe When the probe laser is weak, the continuum resonance is mainly determined by properties of the 
dressing laser. Fluctuations in the dressing laser intensity can disturb the observation of the continuum resonance. 
A change in the dressing laser intensity changes the Stark shifts as well as the peak value and central wavelength of 
the LICS. Because the Stark shifts follow the fluctuations of the dressing laser intensity, the position of the structure 
varies with every laser pulse. The observed LICS is therefore broadened by the fluctuations. In the most extreme 
cases the detected ion signal, for a fixed detuning, fluctuates between the maximum ionization enhancement and 
the ionization suppression. The averaging produced by the fluctuation reduces the LICS contrast between maximum 
enhancement and minimum suppression. 

A weak probe laser does not affect the LICS, and so the probe laser intensity fluctuations have no influence on the 
continuum resonance if there are no dramatic variations in the laser intensity for single pulses. It is not necessary, in 
simulations, to average over intensity fluctuations of a weak probe laser. 

b. Strong probe When the probe laser is strong the continuum is modified by the dressing laser and by the 
probe laser. The influence of the strong probe laser can be observed in the saturation of ionization, i.e. a reduced 
enhancement in the continuum resonance. Although saturation might be reached in the center of the spatial probe 
laser profile, the intensity could still be below saturation in the wings. The maximum enhancement is stronger in the 
wings than in the center of the profile. An integration over the spatial profile of a strong probe laser therefore leads 
to less saturation than expected for constant strong laser intensity. 

Fluctuations in the probe laser intensity determine, as do fluctuations in the dressing laser intensity, additional 
Stark shifts, i.e. variations in linewidth and position of the continuum resonance, as discussed above. 

IV. ANALYTICAL EXPRESSIONS FOR THE LICS PROFILE 

It is always instructive to augment numerical simulation with analyses that can provide explicit formulas for various 
observables. We present in the following sections some examples of such formulas. We consider pulsed excitation, so 
that one cannot always rely on the stationary solutions that are known to give the formulas Eqn. (Q) and Eqn. (|^) 
for the LICS shape. 

A. General theory 

If we neglect the coupling il(t) of Eqn. ([l2|), then the amplitude of C\(t) decays exponentially with time, while its 
phase depends on the Stark shift S\(t). With this in mind we introduce the new variables 
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fti(t) = lnd(t) -i I dt'S x (t'), (19) 

J —OG 

h 2 (t)=C 2 (t)/C 1 (t) (20) 

In terms of the new variables the populations of states 1 and 2 after the interaction with the laser pulses are 

\C x (t)\ 2 = exp[2Re/ii(t)] and \C 2 (t)\ 2 = cxp[2Re/ii(i)] |/i 2 (i)| 2 . (21) 

In the following we assume the probe pulse £ p (t) to be weak, meaning that for the pulse duration r p the direct 
ionization probability 



e 2 



/ di£rg(i) (22) 

«/ oc 



J 

remains small: £ 2 < 1. We therefore seek solutions in the form of a power series in the parameter e. Because for 
e = the populations are |Ci(i)| 2 = 1 and | C2 (t) | 2 = 0, the variables hi(t) vanish for e — and the series expansion 
reads 

k 

We assume that the state 1 is long lived so that one can neglect the spontaneous decay during the interaction with 
the laser pulses (71 = 0) . We assume also that the frequency ui 2 is small enough that the dressing field does not ionize 

the atom in state 1, and hence T\(t) — Y^^ipft)- The equations for h\(t) and fi2(t) derived from Eqn. ( |l2| ) then read 

i 

|M*) = -^ + ^n(t)M*), (23) 



dt 



h 2 (t) = - 



r 2 (t) , 72 r lp (t) 



-iD eff (t))h 2 (t) + 



(iq - 1) 



n(t)(i-h 2 (t) 2 ), 



(24) 



where D e ff{t) = D — S 2 (t) + S±(t) is the two-photon detuning incorporating the Stark-shifted resonance. We require 
solutions that vanish initially: hi(— 00) = h 2 (— 00) = 0. 

The value of h\{i) obtained from Eqn. (23), when integrated over the pulse duration r p , is proportional to e 2 . To 
first order in e Eqn. ( f2~I|) reads 



dt 



h { 2 \t) = - 



^-1-^/(0 ; 



n(t). 



The solution of this equation with the initial condition h 2 (—oo) = is 

dr Q(t — t) exp 



h?(t) = ^ 



1 ,, fr 2 (t') 72 



t-T 



Using Eqns. ( |26|) and ( |23|) we find (t) in second order of perturbation theory to be: 



(25) 



(26) 



Q(t')h 2 1] (t')dt' . 



(27) 



From (26) one can see that h 2 (+00) — and hence, to second order in the small parameter e, the ionization 
probability P is 



P = 1 — cxp(— k), 

where the exponent k is the value of the variable h\(i) in the limit of long times, 



(28) 
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-2Re 



Using Eqn. ( |27[ ) for h^\t) and Eqn. (|26| ) for h^'it) we can write the parameter n as 



k = kq — Re 



-OO OO 



dt J d,Tn(t)n(t-T)exp[-<p(t,T) 

oo 



where Kq is the value of the parameter k in the absence of the dressing laser, 

+oo 



(29) 



(30) 



and ip(t,r) is 



(31) 



Equations ( p^ ) and ( |29| ) include explicitly the time dependence of the effective two-photon detuning, and thus they 
provide the desired analytic formulas for the LICS profile. This expression will be simplified below for some interesting 
limiting cases. 

For small values of k the ionization probability, and hence the observed LICS profile, is proportional to k, 

P Ri K for K« 1. 

In this approximation k provides a theoretical estimate of the observed LICS profile of P; conversely, a parametric fit 
of the detuning dependence of the observed P can be interpreted as the detuning dependence of k. 



B. Short probe pulse 

Maintaining the assumption of a weak probe pulse, Ip ^ <§C il , we now consider laser pulses of different durations, 
and take the probe pulse duration r p to be much shorter than the duration Td of the dressing pulse [r p <C Td). In this 
case the main contribution to the integrals over r and t in Eqn. (^9j) stems from r < t p and \t\ < t p . Hence we can 
write the phase factor of Eqn. (|3l|) in the form 



ip(t, t) 



72 ■ A 

iA 



(32) 



where t m is the time when the intensity of the probe pulse reaches the maximum and A = D e ff(t m ) is the effective 
two-photon detuning at the moment of peak pulse. Using the Fourier expansion of the two-photon Rabi frequency 

n(t) 



f2(w) = — / dui exp (iuit) fl(t) 
2n 



we rewrite Eqn. (E9J) as 



n = kq — 7rRe < (iq — l) 2 / dcu 



[T 2 (tm) + 72]/2 - «A 



(33) 



Equation (J3^) is the sum of two terms. The first one, kq, gives the ionization probability in the absence of the dressing 
laser. The second term involves the convolution of a complex Lorentzian profile with the power spectrum of the two- 
photon Rabi frequency. The shape of this term [i.e. the dependence on the detuning A = D e ff(t m ) at maximum 
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intensity] depends strongly on the relationship between the quantity T2(t m ) + 72 and the bandwidth Alu ~ 1/ T p of 
the two-photon Rabi frequency £l(u>). 

For purposes of illustration we assume a Gaussian shape for the probe and dressing laser pulses 



£i(t) = £ i0 exp 



2r 2 



(34) 



We furthermore assume that only one continuum is involved in the process (this implies that initial and final states 
have angular momentum I = 0). Then Eqn. (33) can be written as 



K 



1 - Re 



(^-l) 2 r 2 



(r d /2)exp [-(cotq) 2 ] 
(T d +72V2- iA + ico 



(35) 



where the time scale r is defined by 



1/t 2 = 1/t 2 + 1/t 2 . 

In Eqn. ( |35| ) r<j is the ionization rate of the state 2 with the dressing field, for t = t m = 0. We shall refer to Eqn. 
( |35| ) as the finite-bandwidth approximation. 

For further simplification we neglect spontaneous emission from state 2 (for very short pulses 72 t p <C 1) and assume 
the dressing pulse to be long, so one can set tq = t p . We obtain 



= 1 - Re 



(jg - l) 2 T p 



dio 



(r rf /2)exp[-(cjr p ) 2 
F d /2 - iA + iu 



(36) 



Eqn. (p6[), taken with Eqn. (Pq), provides an analytic expression for the shape of LICS in the pulsed fields. It is a 
superposition of two profiles: convolution of two symmetric profiles, a Lorentzian and a Gaussian (i.e. a Voigt profile), 
proportional to(l — q 2 ); and a convolution of an asymmetric dispersion-like profile with a Gaussian one, proportional 
to 2q. Figure |l3| shows examples of LICS calculated using Eqn. ( |36|) for different TdT p . We used q = 0.73, appropriate 
for the helium states relevant to our experiment. 

For short pulses (r p < l/r<j) the LICS is broad, with only minor (if any) enhancement and a small dip. For longer 
pulses the LICS is more pronounced. Full suppression is approached for t p > lO/L^. Figure |l4| shows the dependences 
of A m i n r and A max r on the dimcnsionlcss parameter TdT p for q — 0.73. For TdT p < 10 the pulsed nature of the probe 
laser affects the frequencies at which the dip and enhancement of the LICS profile are observed. For a pulsed probe 
laser one can not obtain complete suppression of the ionization probability (see Fig. |l5| ). However, for TdT p — 10 the 
ionization rate of the center of the dip is only 1% of the rate far from resonance. 



C. Corrections to stationary LICS 

It is known [ |39| that for the stationary limit the effective detuning A = A m i n of the dip and of the enhancement, 
A = A max , are 

A m i n — ^"2"' ^max — ~2q' (37) 

The separation w between the dip and the enhancement is proportional to the ionization rate Td, modified by the 
Fano q parameter: 

to = A min - A rnax = ^ (a + iV (38) 

For pulsed excitation the stationary results in Eqns. ( |37j ) must be modified. For TdT p ^> 1 one can find perturbation- 
theory solutions for A m j n and A max which replace Eqn. ( |37| ) with the approximations 

Amin = 9 T ( 1+ (q 2 + l)(rr p ) 2 )' Amax = -^ ( 1+ (g2 + 1)(rdTp)2 )' (39) 
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showing that the width w of the LICS is dependent on the probe-pulse duration r p , 

In this case the values of the parameter n for the effective detuning A = A m i n (n — K m i n ) and for A = A max 
(k = K max ) are determined by the expressions 

Kmm = ? (AD 

2 i 1 \fr ,-r- \2 y > 



k (q 2 + l)(T d T p ) 

— = (<i 2 + i) ( i - rrxw ■ ( 42 ) 

Equations ( 07^ - (^3) show how the transient nature of pulsed LICS affects the profile, raising the minimum value of 
the ionization probability in the dip, lowering the enhancement, and increasing the width of the profile. 



D. The quasistationary case 



The previous analysis shows that the influence of the probe laser pulse duration t p (or the transit time of atoms 
across a probe beam) on the shape of the LICS is important if the dressing-pulse intensity is not very large, T d T p < 10. 
For larger Yd one can neglect the delay in the formation of coherence during the pulses. This means that we can 



(t) — in Eqn. (p5|) and can write, for any r p , 



(iq-l)Q(t) 



T 2 {t)+ l2 -i2D eff {tY 



Then Eqn. (|29|) for the parameter k becomes 

k = k — Re 



)+ l2 -i2D eff {t>) 



(43) 



(44) 



This equation involves a weighted time integration, rather than the frequency integration of Eqn. (|33J). It presents 
the LICS profile as a time-averaged Fano-Cooper or Breit-Wigner profile. The observed profile is affected by the time 
dependence of the width T 2 (t) and the effective detuning D e ff(t) during the interaction. 



E. Stationary case, multiple continua 

Note that even in the case of a time-independent dressing-laser field the LICS profile parameterized as in Eqn. (|^) 
must include the term tTf, when coupling occurs to more than one photoionization continuum. To illustrate this point, 
we assume for simplicity that the probe pulse is much shorter than the dressing pulse, r p r^. We also neglect 
spontaneous emission from state 2, so that 72 = 0. Then Eqn. (Q) can be written in the form of Eqn. (|^) with 

r d = a ^r«(t m ), r = (i-a)£r«(t m ). (45) 
j 3 

The parameter a, describing the effect of multiple continua, is 



jr M *(Erg(t)J (Er2(t) 



< 1- (46) 



The shape of the LICS depends strongly on the relationship between the ionization rates T 3 lk it). If a ^ 1 the enhance- 
ment of the ionization is smaller and the ionization dip does not reach zero as it does for only one continuum, see 



14 



Eqn. (p4[). The minimum and maximum values of n provide information needed to determine the parameters a and 
q, using the formulas 

Kmax/Ko = 1 + aq 2 . 

The parameter a will be unity if there is a single continuum, as occurs for helium when the bound states have orbital 
angular momentum t = 0. (In this case the angular momentum of the continuum states is £ = 1.) Alternatively, 
one has a ~ 1 for a transition within the triplet manifold, 3 Si - 3 Si. In this case three continua are involved in the 
process: 3 Po, 3 Pi, and 3 P2- With the neglect of spin interactions, the radial matrix elements for the bound-continuum 
transitions do not depend on the angular momentum j; with this approximation one finds a — 1 from Eqn. (p6|). In 
all other cases a < 1. 



F. Quasistationary case, single continuum 

Let us consider the case a = 1, so that fl 2 (t) = Tx(t)T2(t) (only one continuum is involved in the process). Then 
Eqn. (^) for the parameter k can be written as 

,«,«,fe|±#,. (48, 
where the variable x(t) is the ratio of a time varying Stark-shifted detuning and a time varying width ^(i). 

Eqn. ([l8|) gives the LICS profile as the integral of a time-dependent single-continuum profile f(x) = (x + q) 2 /(x 2 + 1) 
weighted by the variable T\(t). Consequently the resulting lineshape, Eqn. ( fl8| ) is different from the function f(x): 
the dip does not drop to zero, and the enhancement is smaller. The value of k in the dip is determined by the Stark 
shift and power broadening. It can be estimated for Gaussian pulses and t p <C Td to be 



2 ri 



Kmm l { g + 2[5 2 (0)-S 1 (0)]/r 2 (0)} .„ 
n ~2 q 2 + l rf (50j 

This formula shows the minimum value of the photoionization suppression that can be reached for the quasistatic 
case. 



G. Comparison between theories 



Our primary comparison between experiment and theory was made by means of numerical simulation - solving 
the coupled ordinary differential equations, Eqn. ([ij]). It is instructive to compare the results of those computations 
with several analytic approximations discussed in the previous paragraphs. Figure [l6] presents such a comparison. 
The full-line curve shows the results obtained by integrating Eqn. The dashed line is the finite-bandwidth 

approximation of Eqn. (]35|). The dotted line is the quasistationary approximation of Eqn. (p4|). As can be seen, all 
three curves are very close to each other. For these plots the value of TdT p is 6.4. If this quantity were larger the 
results would approach the dotted line (the quasistationary limit) even closer. 



V. CONCLUSIONS 



We have presented a detailed set of observations of LICS profiles obtained when coupling two states of metastable 
helium to the same continuum by means of pulsed lasers. The profiles exhibit both enhanced and suppressed photoion- 
ization, at appropriate wavelengths, for a narrow LICS in an otherwise flat continuum. Photoionization suppression, 
very pronounced for these observations, shows the presence of population trapping. 
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We have described a simple single-electron model potential whose wavefunctions provide remarkably good values 
for atomic parameters. We have shown, in an appendix, that examples of oscillator strengths computed with this 
potential are very close to values obtained using more elaborate methods of computing wavefunctions. 

We have shown that solutions of the appropriate coupled differential equations, using atomic parameters (Stark 
shifts and Rabi frequencies) computed with the model potential, provide excellent fits to our experimental data. 

We have developed analytic approximations to the photoionization probability produced by pulsed laser pairs, and 
we have shown that these can provide a good approximation to the more elaborate and time-consuming solution of 
differential equations. 
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APPENDIX A: ESTIMATING FANO PARAMETERS 



To extract parameters from resonance profiles of either experimental data or simulations we carried out least-squares 
fits. Less elaborate estimates of the parameters can be obtained by inspection of the profiles. 

Consider the profile function proposed by Fano and Cooper to describe an autoionization resonance, with resonance 
frequency ojr and width T, in the presence of nonrcsonant background. The photoionization cross section is expressible 
in the Fano form 

a(w)=a b {u,)+a a ^^-, (Al) 

where q is the dimensionless Fano q parameter and x is a dimcnsionlcss detuning, 

x = (w - w fl )/(r/2). (A2) 

Neglect of the slow variation of o\, with frequency leads to a profile which has a maximum, a max , at u = uj ma x, 
when x = x max = l/q. It has a minimum, (7 min , at lu — 0J m i n , when x = x min = —q. It follows that a measurement of 
the extrema values of the cross section, a max and o m im an d the associated frequencies u) max and u> m i n , together with 
the asymptotic cross section er(oo) obtained far from resonance, can be used to determine the various parameters of 
the profile: 

(A3) 
(A4) 









Oa = 


tr(oo) 




q = 


&max 


— <j(og) 


a(oo) 





(A5) 

Note that if the resonance formula is applied to a case for which <j\, = (i.e. there arc no additional contributions 
to background apart from the continuum in which the resonance is embedded), and if the LICS profile is normalized 
such that u a = 1, then <T min — and er(oo) = 1, so the expression for q becomes 

q = V^rnax - 1. (A6) 

Once q is known, the width parameter T can be determined from the measured frequencies of the extrema, using 
the expression 

r {^max ^mmjlTT" 2T ^{^max ^min) 7T~~T 2\ ' (A7) 

Z(i + q ) (1 + q ) 

The resonance frequency uir is obtainable from the expression 

= 2^- UJmax + U rnin) - ^ - ■ (A8) 



APPENDIX B: THE MODEL POTENTIAL METHOD 



1. The Radial Equation 

When we introduce a separation of variables 

^(r) = R(r)Y(9, 4>) = r _1 P(r)y(e, <p), 



(Bl) 



into the Schrodingcr equation for a particle of energy E, mass m and charge — e moving in a central potential V(r) 
we obtain the differential equation 



*2 / pa T 2 



P(r)Y(0,<j)) = 0. 



(B2) 



17 



For an attractive Coulomb potential of charge Ze the central potential is 



We take Y(9,(f>) to be an eigenfunction of L 2 , either a single spherical harmonic Yg m (9,<fi) or any combination of 
spherical harmonics with common value of £, so that the differential operator L 2 can be replaced by the eigenvalue 
1(1 +1). It is useful, and customary, to introduce the Bohr radius ao as the unit of length and to express energies in 
multiples of the atomic unit of energy Eau 



4ire Q fr 



Eau = 



47rea m(a ) 2 



(B3) 



The result is the radial equation 



dr 2 



£(£ + !) 



- 2V(r) + 2E 



P{r) = 0. 



(B4) 



This equation has bounded solutions for any positive energy, E > 0, while for negative energies it has bounded 
solutions for a set of discrete values, the bound states of the potential V(r) with angular momentum I. 



2. Effective potentials 

A popular approach to the determination of radial functions is the quantum-defect method originally developed 
by Bates and Damgaard |5Q|. It was used for calculating the ionization rates for helium in |5l|-|53|]. This method 
uses solutions to the radial equation for a Coulomb potential, but forces the bound energies to take specified values. 
Unfortunately this requirement yields functions which diverge at small r, and this poses difficulties in the computation 
of Stark shifts and two-photon Rabi frequencies when the laser frequency is large enough to cause ionization. 

As was shown in |]54| the use of a model potential for the calculation of the atomic parameters such as dynamical 
polarizabilities, Stark shifts, multiphoton ionization rates is convenient. This method is based on the pseudopotential 
theory Ipg-p?! . Because an exact analytical expression for the pseudopotential cannot be derived a number of model 
potential was considered. 

3. Coulomb potentials 

When the effective potential is a Coulomb potential for a charge Ze radial equation is the hydrogenic radial equation 

r d 2 £(£+1) 2Z 



, , , , , 2E 

ar z r z r 



P(r) = 0. (B5) 



As written here, the equation for a Coulomb potential involves three parameters: the nuclear charge Z , the angular 
momentum £, and the energy E. Rather than use the energy E, it proves useful to introduce alternative parameters. 
When E is negative useful real-valued parameters are K and is, 

(hK) 2 ( me 2 \ 2 Z 2 Z 

E = - ' ' = - T\ 7T^ or K= • (B6) 

2m \4:Tre h J 2v 2 a v 

The quantity UK is the rms momentum, while 1/v is proportional to the mean value of Z/r. When E is positive the 
useful real-valued parameters are k and 77, 

(hk) 2 { me 2 \ 2 Z 2 Z 

E =- — — = — I — n or k = . (B7) 

2m \47reo«/ 2r\ l a^rj 

The rms momentum is now hk. 

For a purely Coulomb potential, with nuclear charge Z ', the allowed bound state energies are given by the Bohr 
formula, 
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K = - * (B8) 

and are parameterized by the integer principal quantum number n. When the potential deviates from a purely 
Coulomb potential, then the energies no longer fit this pattern. Often, however, they fit the simple pattern 

Z 2 
2(n - d e y 

where 8i is the quantum defect. There are several ways of constructing wavefunctions associated with such energies. 



4. The Radial Functions 



Consider a generalization of the negative-energy hydrogenic radial equation 

2" 



d 2 A(A + 1) 2Z (Z 
dr 2 r 2 r 



P(r) = 



with the energy (in atomic units) defined as 



E 



2v 2 ' 



This radial equation has negative-energy solutions of the form 



P(r) = Wexp - 



Zr 



„A+1 



2Zr 

1 Fi(-i/ + A + l;2A + 2; ) 



where i-Fi(a; c; x) is a confluent hypergeometric function (Kummer function), 

ax a{a + 1) x 2 
1 F 1 (a; C ;x) = l + -- + ^ TTj - + ... 

T(a + n) T(c) x n 
T{a) T(c + n)T(n + l) 



n=0 



and where Af is a normalization constant, expressible variously as 

TV 



2Z 
v 



2y/Z jdE T(v + \+l) 
r(2£ + 2)Vd^y I>-A) 

2Z\ e 2Z 3 / 2 1 jr(n r + 2A + 2) 
v 2 T(2£ + 2)\ vT(n r + l) 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



where n r — v — A — 1 is the radial quantum number. The radial functions can also be written in terms of Laguerre 
polynomials 



/v> = (- A ' n "' +1 ' 



x exp 



/' / v 2 y T(n r + 2A + 2) 

Zr \ r A+l L (2A+l) f^ r ^ 
V J Ur \ V 

The energy-normalized radial wavefunctions for positive energy E — k 2 /2 can be written as [ ]60[ 

|r(A + l + i?7)|exp (7777/2) 



(B15) 



P E e(r) 



(2kr) 



A+l 



T(2A + 2) V2t7I 
x exp (-ifcr) r A+1 x Fi (A + 1 + iry; 2A + 2; 2z£;r) . 

where A: = V2E is the mean momentum and 77 = Z/fc is the Coulomb parameter. 



(B16) 
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5. Examples 



c. Hydrogenic Ions When A is an integer, A = £, as is the case for hydrogen, then v must be an integer, the 
principal quantum number n, in order that the hypergeometric function remain bounded: 

if \ = £ then v = n. (B17) 

This condition fixes the energy E by the Bohr formula 

K = ~- (B18) 

d. Screened Hydrogenic Functions One method to obtain a wavefunction for one or two low-lying bound states is 
to introduce an effective nuclear charge Z* through the definition 

Z* nl = ^-2E ne . (B19) 

This approach maintains the quantum numbers n and £, and deals with well behaved bound functions. 

e. The Quantum Defect Method One approach to obtaining radial functions with the desired energies is to simply 
let v be defined to produce the desired energies, 

u = v n =n-6 t , (B20) 

while retaining the integer values for £. However, this approach, commonly termed the quantum defect method, does 
not produce radial functions which remain bounded for all distances. In order that the hypergeometric function be a 
polynomial, it is necessary that the first argument be an integer, the radial quantum number n r , 

n r = A + 1 - v. (B21) 

/. The Effective Potential However, we can choose A to be a parameter fixed by £ and a given energy: we write 
the (integer) radial quantum number as 

n r = -n + £+l = -u + X+l. (B22) 
This means that the energy is given by the quantum defect formula 

^ = - 2 { n + Z l-£f (B23) 

and, in turn, that the parameter v is 



v <* = y-^= n + X - i - ( B24 ) 

One can use the lowest energy value E\ (i.e. the ionization energy) to define the quantum defect, so that one has, for 
specified £, the effective orbital angular momentum 



X = l-l + J-£r. (B25) 



Then the nth radial function for specified £ is 



7r\ 9Zr 
P{r) = Mr x+1 exp iF^-n + £ + 1; 2A + 2; ). (B26) 
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6. Dipole Transition Moments 



The radial matrix elements of the dipole transition moment for transition between states with the radial quantum 
numbers n a and n b and with effective angular momentum £ a and l b respectively is 



R a>b = I P a (r)P b (r)rdr. 
lo 



(B27) 



Many authors have presented expressions for these matrix elements when the radial functions are those of hydrogen-like 
atoms f6l|-|72]| . We have found that a convenient expression is 

1 ( 2v b \ Xa+2 ( 2v a \ Xb+2 r(A a + A fc + 4) 



±Z\v a +Vb) \v a + v b J r(2A a + 2)r(2A 6 + 2) 

r(i/ a + A + i)r(i/ b + A 6 + i) 



-,1/2 



T(n a + l)T(n b + 1) 



X F 2 [X a + X b + 4; -n a , -n b ; 2A a + 2, 2A h + 2; 



2v b 2v a 



Va + V b V a + V b 



(B28) 



where F 2 is the generalized hypergeometric function of two arguments |73| ] : 

T{a + n + to) r(/3i + n) T{(3 2 + m) 



E 



n.m—O 



T(a) r(/3i) r(/3 2 ) 



r(7i) r( 72 ) 



V" 



r( 7 i + n) T( 72 + n) r(n + 1) r(m + 1) ' 



(B29) 



When /3i = — n a and /3 2 = —n b are negative integers the series terminates and F 2 is then a polynomial in x and y. 
For transitions from a bound state a to the continuum state b one can write 



Ra,b= I P a (r)P b (r)r dr 
lo 



|r(A b + l + zZ/fc fc )|r(A a + A h +4) 
(2fc) 3 r(2A a + 2)r(2A 6 + 2) 



k b T {v a + A a + 1) 



-1 1/2 



27rZr(n a + 1) 



2kv a \ Xb+3 ( 2Z \ Aa+1 2 ( vZ 



Z + ik b v a 



Z + ik b v„ 



— ex P \ .,, 



Z" 2Z 2ik b v n 
x F 2 A a + A h + 4; -n», A fc + 1 + 2 — ; 2A a + 2, 2A b + 2; — — , b 

Kfc Z + Z + Z/Cfel'a 



(B30) 



In this case the function F 2 is the sum of finite number of terms that are proportional to the hypergeometric function 
of one argument 



2ik b v a 
Z + ik b v a 



(B31) 



Using the expression ( B29| ) and (B30) one can compute, with sufficiently high accuracy, values of the ionization 
rates, Stark shifts and two-photon Rabi frequencies. 
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7. Needed parameters 



For Raman transitions between S states (involving only one P continuum) the ionization rates of Eqn. ( |l5|) , Stark 
shifts of Eqn. ( |l6| ) and two- photon Rabi frequency of Eqn. (18) can be expressed in simpler form fl^| : 



T nj = (4^ 2 /3) alj \f{n,io 3 + E n \ 2 , 



(B32) 



S nj = (4tt/3) alj \9(n, m)f 



E„ 



E„ - uj. 



(47r/3)a/j dE\f(n,E)\' 



1 



-E + E n + uji -E + E n 
where a is the Sommerfeld fine-structure constant. The Fano q parameter is 

1 1 



1 

Q = ~ 

7T 



J2g{l,m)g(2, 



m ) 



E„ 



Ei — ujp E m — E\ + u>d 



(B33) 



dEf(l,E)f(l,E) 



1 



1 



E — Ei — Lu p E — Ei — ujd 



[f(l,u;i+Ei)f(2,u J2 + E 2 )} 



(B34) 



Here g(n,m) = R n o m i, f(p,E) = R n o ei are the radial matrix elements determined by the expressions (B28) and 

(B- 

One can see that for the calculation of the Stark shifts and the Fano q parameter we must find the infinite sums in 
(B33) and (B34). In the present calculations we usually restrict the sum to the first 200 terms. This gives an accuracy 
of the infinite sum that is in any case much greater than the accuracy of the model potential method. The integral 
over the continuum spectrum after the transformation is [B4j 



/(*) 



dx 



f(x) - f(a) 



dx 



fix) 



dx 



(B35) 



was computed by Gauss quadrature. 
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APPENDIX C: NUMERICAL VALIDATION OF MODEL POTENTIAL 



1. The Wavefunctions 



Theoretical prediction of a LICS structure rec 
Eqn. (fL6|) and two-photon Rabi frequency Eqn 



uires computation of the ionization rates Eqn. (|l5| ) , Stark shifts 
18|) using realistic wavefunctions. Because these various parameters 
all require summations over many discrete states and integration over many continuum states, it is desirable to have 
wavefunctions which require as little computational effort as possible, while still providing accurate depiction of the 
essential properties of the active electron. We have found that a simple single-electron model potential satisfies both 
these requirements. Appendix [b| describes the mathematics of this method for computing matrix elements. As noted 
there, a single empirically determined parameter Xg fixes the potential, and hence all the wavefunctions, for each 
orbital angular momentum I. 



2. Calculations of the parameters of He 

Because the singlet and triplet systems of He interacts only weakly, we consider each system to be independent. 
For each system we calculate \g using the experimental term values from the NBS tables [Q. Some results of these 
calculations are shown in table IV and in table As one can see from these tables, the value of A; depends weakly 



on the energy. The radial transition moments are also slowly varying functions of energy. However, the Stark shifts 
and two-photon Rabi frequency are very sensitive to the energy of intermediate states: when an accidental resonance 
occurs then the use of an approximate energy value determined by v n can lead to large inaccuracy. 

TABLE IV. The values of the energies, parameteres Xi and v„ for singlet He 





E, cm -1 








E, cm" 1 


X e 


Vnl 




E, cm -1 


A; 


Vnl | 


2 l S 


-32032.3 


-.149 


1.851 


2 X P 


-27175.8 


1 


010 


2.010 












-13445.9 


-.143 


2.857 


3 X P 


-12101.5 


1 


011 


3.011 


to 


-12205.7 


1.998 


2.998 


4 1 5' 


-7370.5 


-.141 


3.859 


4 X P 


-6818.1 


1 


012 


4.012 


4 X D 


-6864.2 


1.998 


3.998 




-4647.2 


-.141 


4.859 


5 X P 


-4368.3 


1 


012 


5.012 


5'D 


-4392.4 


1.998 


4.998 


6*5 


-3195.8 


-.140 


5.860 


6 1 P 


-3035.8 


1 


012 


6.012 


Q l D 


-3049.9 


1.998 


5.998 


7 1 S 


-2331.8 


-.140 


6.860 


7 1 P 


-2231.6 


1 


013 


7.013 


7 1 D 


-2240.6 


1.998 


6.998 


8^ 


-1777.0 


-.139 


7.861 


8 1 P 


-1709.4 


1 


012 


8.012 


d}D 


-1715.2 


1.998 


8.998 


9 X S 


-1397.9 


-.140 


8.860 


9 X P 


-1351.1 


1 


012 


9.012 


9 1 !) 


-1355.5 


1.998 


9.998 


lO 1 ^ 


-1128.6 


-.140 


9.860 


10 1 P 


-1094.8 


1 


013 


10.013 


lO 1 ^ 


-1097.9 


1.998 


9.998 










15 X P 


-486.9 


1 


013 


15.013 


















20 1 P 


-274.0 


1 


013 


20.013 










cont. 




-.139 




cont. 




1 


013 




cont. 




1.998 
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TABLE V. The values of the energies, parameteres A^ and v n for triplet He 





E, cm' 1 




v n l 




E, cm" 1 


A, 


V n l 




E, cm- 1 


Af 


Vnl | 


2 3 S' 


-38454.7 


-.311 


1.689 


2 3 P 


-29223.8 


.938 


1.938 










3 3 S 


-15073.9 


-.302 


2.698 


3 3 P 


-12746.1 


.934 


2.934 


3 3 D 


-12209.1 


1.998 


2.998 


4 3 S" 


-8012.5 


-.299 


3.701 


4 3 P 


-7093.5 


.933 


3.933 


4 3 D 


-6866.2 


1.998 


3.998 


5 3 S 


-4963.6 


-.298 


4.702 


5 3 P 


-4509.9 


.933 


4.933 


5 3 D 


-4393.5 


1.998 


4.998 


6 3 S 


-3374.5 


-.297 


5.703 


6 3 P 


-3117.8 


.933 


5.933 


6 3 D 


-3050.6 


1.998 


5.998 


7 3 S 


-2442.4 


-.297 


6.703 


7 3 P 


-2283.3 


.933 


6.933 


7 3 D 


-2241.0 


1.998 


6.998 


8 3 S 


-1849.2 


-.297 


7.703 


8 3 P 


-1743.9 


.933 


7.933 


8 3 D 


-1715.6 


1.998 


8.998 


9 3 S 


-1448.6 


-.296 


8.704 


9 3 P 


-1375.3 


.933 


8.933 


9 3 D 


-1355.4 


1.998 


9.998 


W 3 S 


-1165.2 


-.296 


9.704 


10 3 P 


-1112.4 


.932 


9.932 


10 3 D 


-1097.7 


1.999 


9.999 


15 3 S 


-508.4 


-.308 


14.692 


15 3 P 


-491.9 


.936 


14.936 


















20 3 P 


-275.9 


942 


19.942 










cont. 




-.300 




cont. 




0.940 




cont. 




1.998 





3. Accuracy 



To demonstrate the accuracy of the model-potential method we display in tables VI and VII the computed oscillator 
strength for S-P transitions in singlet and triplet He. There is exceptionally good agreement between our values and 
those obtained with more sophisticated methods. Except for transitions from the ground state, the differences between 
our values and the more accurate ones are always less than 10%. The lack of accuracy for transitions involving the 
ground state is expected because here the influence of polarization and exchange effects are greatest. 

To test the accuracy of the continuum wavcfunctions we have calculated the cross sections for photoionizations 
from some excited levels. Comparisons of the cross sections calculated with the help of model potential with results of 
other author are illustrated in Fig. [l7|, for the 2s 1 So level. One can see that our values, obtained with a very simple 
approach, are in very good agreement with previous ones obtained by more sophisticated methods. The accuracy of 
our one-photon bound-bound and bound-free transitions validates the use of model-potential wavefunctions for the 
calculation of Stark shifts and two-photon Rabi frequency for helium. 
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TABLE VI. Oscillator strength for transitions between 1 S and 1 P-states of He 







2"P 


3*P 


4"P 


5 X P 


vs 


a 


0.362 


0.082 


0.031 


0.015 




b 


0.276 


0.073 


0.030 


0.015 




c 


0.276 


0.073 


0.030 


0.015 


2 1 S 


a 


0.375 


0.165 


0.053 


0.024 




b 


0.376 


0.151 


0.049 


0.022 




c 


0.376 


0.151 


0.049 


0.020 


4 X S 


a 


0.024 


0.291 


0.859 


0.157 




b 


0.026 


0.308 


0.858 


0.146 




c 


0.026 


0.306 


0.855 


0.150 


5 X S 


a 


0.009 


0.053 


0.460 


1.09 




b 


0.010 


0.056 


0.476 


1.08 




c 


0.010 


0.055 


0.470 


1.10 



a The model potential method. b The results of Chen }48|. c The results of SchifT et al 



TABLE VII. Oscillator strength for transitions between 3 S and 3 P-states of He 







2 a P 


3 3 P 


4 3 P 


5 3 P 


1 3 S 


a 


0.563 


0.060 


0.024 


0.012 




b 


0.539 


0.064 


0.026 


0.012 




c 


0.539 


0.064 


0.026 


0.013 


2 3 S 


a 


0.210 


0.908 


0.051 


0.023 




b 


0.209 


0.891 


0.050 


0.023 




c 


0.209 


0.891 


0.050 


0.023 


4 3 S 


a 


0.033 


0.438 


1.225 


0.044 




b 


0.032 


0.436 


1.215 


0.044 




c 


0.032 


0.436 


1.215 


0.044 


5 3 S 


a 


0.012 


0.068 


0.669 


1.540 




b 


0.011 


0.068 


0.668 


1.531 




c 


0.011 


0.068 


0.670 


1.530 



a The model potential method. b The results of Chen c The results of Schiff et al 
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FIG. 1. Laser induced couplings in metastable helium. The probe laser at 294 nm and the dressing laser at 1064 nm couple 
the metastable singlet state 2s 1 So and the initially unpopulated excited state 4s So through the ionization continuum. For 
strong laser intensities two-photon ionization (IR and UV) of the metastable triplet state 2s 3 Si forms an additional background 
for the ion signal. 

FIG. 2. Layout of experiment, showing atomic beam, discharge, skimmer and lasers. The helium atoms are excited in a 
pulsed gas discharge. The atomic beam then interacts with the lasers in the center of a capacitor. The ions produced there are 
accelerated and detected on a microsphere plate. 

FIG. 3. REMPI transitions for detecting population in the metastable states of helium. The singlet state 2s 1 5*0 is coupled 
with pulsed radiation at 397 nm through the excited state 4p Pi to the ionization continuum. The metastable triplet state 
2s 3 5*i is ionized with laser light at 389 nm through the states 3p 3 Po,i,2- 

FIG. 4. Ion signal resulting from the ionization of the metastable states after excitation of the helium atoms in a pulsed 
gas discharge. The populations are probed by resonant enhanced multiphoton ionization (REMPI) with a pulsed dye laser 
(bandwidth 6 GHz). The probing transitions are strongly saturation-broadened. In the gas discharge about 90 % of the excited 
helium atoms are found in the metastable triplet state 2s 3 Si, and only 10 % are in the singlet state 2s 1 So. 

FIG. 5. Difference between ion signals with and without IR laser. The probe laser frequency tuned far away from the 
two-photon resonance where LICS is observed. The pump laser intensity is varied, while the dressing laser intensity is kept 
constant at 1^ — 226 MW/cm 2 . The measured data are normalized to the total ionization signal from singlet- and triplet 
state. Good agreement is found between the experimental values and the calculated results (dashed line). 

FIG. 6. LICS for coincident (left side) and delayed (right side) laser pulses, when the dressing laser intensity is varied. The 
experimental data are compared to results from numerical simulations (dotted lines), including laser intensity fluctuations and 
integration over the spatial profile of the pump laser. Enhancement as well as strong suppression of ionization are observed. 

FIG. 7. Observed and calculated ionization probability of state 2s 1 So, as the pump laser intensity is increased. The dressing 
laser is switched off. 

FIG. 8. LICS profile for weak and strong pump laser. The dressing laser intensity is 7^ = 44 MW/cm 2 (Fig. 8 a). The 
probe laser intensity is 7p ' = 50 MW/cm 2 (Fig. 8 b). Both laser pulses are coincident. With increased pump laser intensity, 
the maximum enhancement as well as the suppression of ionization is reduced. 

FIG. 9. An example of fitting of the experimental data using Eqn. (^). Squares are experimental data, dotted line is 
the result of fitting, circles are the difference between experimental data and the fit (At — ns, Ip°^= 4 MW/cm 2 , I d °^ = 75 
MW/cm 2 ). 

FIG. 10. Width parameters vs. intensity of the dressing pulse (see Eqn. (^|) and Table |): (1) is the variation of the overall 
width parameter T, (2) is the empirically determined F<j, (3) is the calculated photoionization rate Td- 



FIG. 11. LICS profile for two approximations: the solid line shows the simulation results for a single pair of pulses, the 
dotted line shows the profile assuming steady-state conditions (q = 0.73, 7^ 0) = 2.5 MW/cm 2 , I d 0) = 38 MW/cm 2 ). 

FIG. 12. Dressing pulse and delayed probe pulse (At = -5 ns, 7^ 0) = 4MW/cm 2 , I d 0> = 75 MW/cm 2 ). 
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FIG. 13. Plots of k/ko versus D/Td (detuning in units of the width) for q — 0.73 and various values of TdT p . For given Td 
this product is proportional to the length of the probe laser pulse. 

FIG. 14. Probe laser frequencies at which occurs the minimum (1) and the maximum (4) in the LICS profile of ionization 
P versus TdT p . The effect of the duration of the probe-laser pulse (with long dressing pulse) is clearly revealed. Steady-state 
values (2 and 3) are shown by the dashed lines. They approach the former for large YdT v . 



FIG. 15. Values of maximum (1) and minimum (2) in the LICS profile of ionization P versus YdT p showing effect of pulse 
duration of probe pulse. 



FIG. 16. LICS profile in various approximations. The solid line shows the result obtained by numerically solving Eqn. (y_2J). 
The dashed line is the finite-bandwidth theory, Eqn (^). The dotted line is the quasistationary theory, Eqn (jiip. Here the 
dressing laser intensity is I d = 38 MW/cm 2 , the probe laser intensity is I p = 2.5 MW/cm 2 , q — 0.73, r v — 2.3 ns, Td = 5.1 
ns, the laser pulses are coincident. 

FIG. 17. Photoionization cross sections for the 2 1 S state of He: 1 - this work, 2 - @, 3 - f^, 4 - @ ( lMb=10~ 18 cm 2 ). 
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